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ABSTRACT 
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Goda's  (1975)  model,  describing  wave  transformation  from 
deep  water  to  across  the  surf  zone,  is  compared  with  a  large 
amount  of  wave  data  obtained  from  experiments  conducted  at 
Torrey  Pines  Beach,  San  Diego,  California.  Goda's  model 
simulates  wave  breaking  by  truncating  the  Rayleigh  distribu¬ 
tion  in  order  to  estimate  the  wave  height  distributions  across 
the  surf  zone;  wave  heights  are  shoaled  by  applying  nonlinear 
theory.  Comparisons  between  the  empirical  distributions  and 
theoretical  distributions,  and  between  measured  and  theoreti¬ 
cal  rms  wave  heights,  are  made.  It  is  found  that  Goda's 
model  over-predicts  the  tails  and  under-predicts  the  peaks 
of  the  empirical  distributions,  and  that  the  calculated  rms 
wave  heights  are  too  large  compared  with  measured  values. 

The  range  of  breaking,  and  the  coefficients  used  in  the 
breaking  criteria  by  Goda,  are  modified  in  order  to  obtain 
a  model  which  better  fits  the  distribution  of  observed 
heights,  and  which  matches  the  model  and  observed  rms  wave 
heights.  The  results  are  quite  good,  with  error  envelope 
for  predicted  rms  wave  heights  less  than  20%.  Linear  shoal¬ 
ing  theory  is  applied  to  the  model  and  found  to  be  as  good 
as  applying  nonlinear  theory. 
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I .  INTRODUCTION 


The  evaluation  of  an  irregular  group  of  shoaling  waves 
as  they  approach  and  pass  through  the  breaker  zone  is  a  com¬ 
plex  process  which  requires  special  measurements  and  analy¬ 
sis  considerations.  The  usual  approach  to  shallow  water 
wave  transformations  is  to  predict,  from  a  single  "represen¬ 
tative"  set  of  deep  water  parameters ,  the  wave  height,  the 
wavelength  and  the  frequency  at  specific  shallow  water  depths, 
using  linear  shoaling  theory.  The  primary  objection  to  this 
approach  is  that  a  single  set  of  deep  water  wave  parameters 
does  not  realistically  represent  the  distributional  charac¬ 
teristics  of  naturally  occurring  sea  surface  waves.  A 
secondary  objection  arises  from  the  use  of  linear  transfor¬ 
mations  which  become  inadequate  when  applied  through  the 
surf  zone  (Wood,  1974) . 

Wave  heights  in  deep  water,  having  Gaussian  surface  ele¬ 
vations,  are  described  by  the  Rayleigh  distribution  (Longuet- 
Higgins,  1952).  Waves  propagating  toward  shore  can  increase 
in  height  due  to  shoaling  effects,  refraction  and  wave  inter¬ 
actions  ,  and  eventually  reach  a  depth  where  they  start  break¬ 
ing.  The  energy  dissipation  due  to  breaking  has  been  simu¬ 
lated  (Goda,  1975)  by  truncating  the  tail  of  the  Rayleigh 
distribution. 

Experiments  were  conducted  at  Torrey  Pines  Beach,  San 
Diego,  California,  during  November  1978.  Sea  surface  eleva¬ 
tions,  pressures  and  velocities  were  measured  at  closely 
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spaced  locations  in  a  line  extending  from  10m  depth  to 
inside  the  surf  zone.  This  thesis  applies  Goda's  model  to 
the  measurements  in  order  to  examine  the  shoaling  and  trans¬ 
formation  of  wave  heights  and  their  probability  density  func¬ 
tions  (pdf's)  from  deep  water  to  breaking  and  across  the 
surf  zone  to  the  shoreline. 
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II.  THEORETICAL  BACKGROUND 


A.  RAYLEIGH  DISTRIBUTION 

The  Rayleigh  distribution  was  shown  theoretically  by 
Longuet-Higgins  (1952)  to  apply  to  deep  water  wave  heights 
on  the  assumption  that  the  sea  waves  are  a  narrow-banded 
Gaussian  process.  Barber  (1950)  had  earlier  presented  empiri¬ 
cal  evidence  that  the  Rayleigh  distribution  agreed  with  the 
measured  distribution  of  waves .  On  the  assumption  that  the 
wave  height  is  twice  the  wave  amplitude,  the  wave  height 
probability  density  is  then  represented  by 

pun  -  2H/H^s  expl-H2/^)  (1) 

where  H _  is  the  rms  wave  height. 

rms 

Using  pressure  records  in  the  Gulf  of  Mexico,  Longuet- 
Higgins  (1975)  observed  that  the  Rayleigh  distribution  fits 
the  observed  distribution  reasonably  well  in  "fairly  deep 
water" .  He  found  that  there  is  a  slight  excess  of  waves 
with  heights  near  the  middle  of  the  range  and  a  deficit  at 
the  two  extremes.  Since  much  of  the  high-frequency  portions 
of  the  wave  records  were  filtered  out  by  the  pressure  trans¬ 
ducer,  Longuet-Higgins  (1975)  suggested  that  the  narrow  band 
approximation  may  not  be  as  applicable  for  the  unfiltered 
records.  In  shallow  water  with  much  steeper  waves,  the 
Rayleigh  distribution  can  again  be  expected  to  be  less 
applicable  due  to  the  non-linearities. 


Since  the  Rayleigh  distribution  theoretically  did  not 
apply  to  broadband  wave  spectra,  Goda  (1970)  numerically 
simulated  wave  profiles,  where  the  amplitudes  were  specified 
by  various  theoretical  spectra  of  varying  bandwidth  and  the 
phase  was  random.  He  then  examined  the  simulated  records  for 
surface  elevations,  crest-to-trough  wave  heights  and  zero- 
up-crossing  wave  heights.  He  found  that,  using  the  zero-up- 
crossing  determination  of  wave  heights,  the  Rayleigh  distri¬ 
bution  is  a  good  approximation  irrespective  of  the  spectral 
bandwidth.  Tayfun  (1977) ,  in  studying  the  transformation  of 
deep  water  waves  to  shallow  water  waves,  showed  that  the 
Rayleigh  distribution  for  wave  amplitude  was  generally 
applicable  to  all  bandwidths . 

The  Rayleigh  distribution  is  applied  correctly  only  to 
low  waves  in  deep  water  (Longuet-Higgins,  1975) ,  since  it  is 
assumed  that  the  contributions  from  different  parts  of  the 
generating  area  are  linearly  superposable .  Under  this  assumpt 
tion,  the  distribution  clearly  should  not  hold  for  waves 
approaching  maximum  height,  i.e.,  close  to  breaking,  as  in 
the  surf  region  or  even  in  the  open  sea  with  whitecaps .  It 
has  been  found  by  several  authors  (Chakrabarti  and  Cooley, 
1977;  Forristal,  1978)  that  the  theoretical  Rayleigh  dis¬ 
tribution  over-predicts  the  maximum  wave  in  the  tail  compared 
with  large  wave  observations.  Forristal  (1978)  attributed 
the  differences  to  the  non-linear,  non-Gauss ian  and  skewed 
nature  of  the  free  surface. 


Tayfun  (1980)  examined  non-linear  effects  by  consider¬ 
ing  an  amplitude-modulated  Stokian  wave  process  with  the 
restriction  that  the  underlaying  first  order  spectrum  is 
narrow  band.  The  surface  displacements  were  found  to  be 
non-Gaussian  and  skewed,  and  wave  heights  distributed  ac¬ 
cording  to  the  Rayleigh  probability  law,  particularly  for  low 
and  medium  wave  height  ranges.  On  the  basis  of  the  results 
obtained,  Tayfun  (1980)  concludes  that  the  non-Gaussian 
characteristics  of  the  free  surface  do  not  directly  result 
in  reducing  maximum  wave  heights  in  a  manner  consistent  with 
field  observations  and  that  a  more  plausible  mechanism  is 
wave  breaking,  which  is  a  non-linear  effect  not  directly 
accounted  for  in  the  analytical  wave  models  currently  available. 

Longuet-Higgins  (1980)  analyzed  the  effects  of  non¬ 
linearity  and  finite  bandwidth  on  the  distribution  of  wave 
heights  to  explain  the  differences  with  observations  found 
by  Forristal.  He  found  that  the  reason  for  the  discrepancy 
could  be  accounted  for  by  the  presence  of  free  background 
"noise"  in  the  spectrum,  outside  the  dominant  peak,  and  that 
it  was  not  due  to  a  finite-amplitude  effect.  Longuet-Higgins 
concludes  that  the  distribution  of  wave  heights  even  in  a 
storm  is  well  described  by  the  Rayleigh  distribution,  pro¬ 
vided  the  rms  amplitude,  a,  is  estimated  from  the  original 

—  1/2 

record  and  not  from  the  frequency  spectrum  as  a  =  (2mQ)  '  . 

The  effect  of  finite  bandwidth  is  estimated  from  a  model 
assuming  low  background  noise  linearly  superposed  on  a  very 
narrow  (delta  function)  spectrum.  For  narrow  bandwidths. 
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he  obtains  the  formula 


a2/2m0  =  1  -  0.734  v2  (2) 

where  v  is  the  rms  spread  of  the  noise  about  the  mean  fre- 

2 

quency.  Values  of  v  corresponding  to  Pierson-Moskowitz 
(broad-band)  spectra  also  give  results  in  close  agreement 
with  observation.  Therefore,  the  Rayleigh  distributions 
calculated  in  this  paper  are  parameterized  using  the  rms 
wave  height. 

B.  TRUNCATED  PROBABILITY  DISTRIBUTIONS 

In  concept,  waves  are  described  by  the  joint  distribu¬ 
tion  of  height,  period  (or  equivalent  wavelength)  and  direc¬ 
tion.  To  simplify  the  analysis,  all  authors  assume  a  very 
narrow  band  frequency  spectrum  and  a  narrow  directional  spec¬ 
trum,  so  that  all  the  wave  heights  of  the  distribution  are 
associated  with  a  single  mean  frequency  and  mean  direction. 
Therefore,  starting  in  deep  water,  the  waves  are  described 
by  the  unaltered  single-parameter  Rayleigh  distribution, 
with  the  implied  assumptions.  The  deep  water  wave  heights 
are  transformed  into  shallow  water  waves  using  shoaling 
theory  in  which  frictional  dissipation  is  neglected.  Even¬ 
tually  the  waves  reach  such  shallow  water  that  they  start 
to  break,  with  the  largest  waves  breaking  furthest  offshore 
first.  Wave  breaking  is  simulated  by  truncating  the  tail 
of  the  Rayleigh  distribution. 


1.  Collins  Distribution 


Collins  (1970)  was  the  first  to  apply  the  technique 
of  a  truncated  distribution  to  describe  the  effects  of  wave 
breaking,  using  a  sharp  cut-off  with  all  broken  waves  equal 
to  which  results  in  a  delta  function  at  H^.  Collins  does 
not  give  an  explicit  formula  for  his  distribution,  but  it 
would  be  the  same  as  the  described  by  Battjes  (1974)  (see 
below) .  He  used  linear  shoaling  theory  and  the  breaking 
criterion  after  Le  Mehaute  and  Koh  (1967 

Hb/HQ  =  0.76  tan1/7$(H0/L0) ~1/4  (3) 

where  tan  3  is  the  bottom  slope  and  Hq  and  LQ  are  the  deep 
water  wave  height  and  length,  respectively.  For  waves  break¬ 
ing  at  an  angle  a,  the  bottom  slope  is  actually  tan  6  cos  ab, 

1/2 

and  Hq  should  be  replaced  by  Hq  cos  '  <Xq  .  The  various  dis¬ 
tributions,  and  how  they  are  truncated  are  showed  schemati¬ 
cally  in  Table  I . 

2.  Battjes  Distribution 

Battjes  (1974,  1978)  again  used  a  sharp  cut-off  of 
waves  and  applied  the  breaking  criterion  based  on  Miche's 
formula  for  the  maximum  height  of  periodic  waves  of  constant 
form. 


1  0.88/k  tanh  (y/0.88  kh)  (4) 

where  y  is  an  adjustable  coefficient.  In  shallow  water  (4) 
reduces  to 
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Hb  =  Y  h  .  (5) 

He  uses  linear  theory  to  shoal  the  waves.  The  probability 
density  for  breaking  wave  heights  is  given  by: 

p (H)  =  H/2H2  exp [-1/2 (H2/H2) ] ,  for  0  <  H  <  Hb  (6) 

p  (H)  =  exp[-l/2(H2/H2)  ]6(H  -Hb)  ,  for  H>Hb  (7) 

A 

where  H  is  the  wave  height  parameterizing  the  truncated 
Rayleigh  distribution  by  Battjes.  All  waves  that  have 
broken,  or  are  breaking,  assume  the  height  prescribed  by 
(7) ;  this  results  in  a  delta  function  at  the  truncation 
height,  Hb,  of  the  distribution  (see  Table  I) . 

3.  Kuo  and  Kuo  Distribution 

Kuo  and  Kuo  (1974)  investigated  the  effect  of  break¬ 
ing  on  wave  statistics  using  a  conditional  Rayleigh  distri¬ 
bution  sharply  truncated,  specified  by  the  breaking  wave 
height  simply  proportional  to  local  water  depth,  equation 
(5) .  The  conditional  probability  density  function  of  wave 
heights,  pb(H),  is  calculated  using  the  following  equation, 

p  (H)  =  p  (H/0  <  H  <  Hb)  =  -  ,  for  0  <  H  <  Hb 

D  b 

/  p(H)dH 

0 

(8) 

=  0  ,  for  H  >  Hb  . 
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Describing  the  conditional  probability  in  this  manner  re¬ 
sults  in  the  proportional  redistribution  of  probability 
density  associated  with  the  broken  or  breaking  waves  over  the 
range  of  H.  This  removes  the  delta  function  at  the  breaking 
wave  height,  H^,  previously  described  by  Collins  and  Battjes. 
Table  I  shows  the  original  Rayleigh  distribution  with  dotted 
lines,  and  the  modified  Rayleigh  distributions  in  heavy  lines 
after  applying  the  cut-off  to  the  tails  using  the  breaking 
criterion.  The  distribution  by  Kuo  and  Kuo  is  more  realis¬ 
tic  but  still  results  in  a  sharp  cut-off  of  the  distribution 
at  the  breaking  heights . 

5 .  Goda  Distribution 

Goda  (1975)  derived  a  more  realistically  truncated 
distribution,  qualitatively  anyway,  by  requiring  a  gradual 
cut-off  of  the  distribution.  He  uses  a  shoaled  Rayleigh 
distribution  to  describe  the  unbroken  wave  heights  at  shore¬ 
ward  locations  prior  to  applying  his  cut-off  which  is  given 
by 


p  (H)  =  4H/K2H2  exp(-2H2/K2H2)  (9) 

o  so  so 

where  K  is  the  shoaling  coefficient.  Goda  (1975)  calculated 
s 

the  wave  shoaling  using  the  nonlinear  theory  of  Shuto  (1974)  : 
which  dictates  the  following: 


2  2 

0  <  gHT  /h  <_  30:  Small  Amplitude  Theory 


(10) 


30  <  gHT2/h2  <_  50:  H  h2/7  = 


constant 


(11) 


1 

50  <  gHT2/h2  <  •:  Hh5/2 [  v/gHT2/h2  -  2/3]  =  constant.  (12) 

Goda  assumes  that  wave  breaking  occurs  in  a  range  of  wave 

heights  between  H0  and  H.  ,  with  varying  probability.  The  1 

^  x  j 

probability  density  function  of  unbroken  waves  only  is 
expressed  as: 

Pr(X)  =  PQ(X)  ;  for  X  <  X2 

x  -  x 

Pr(X)  =  Pom  -x_^po(X1);  forX2<X<Xl 

Pr(X)  =  0  ;  for  Xx  <  X 

where,  normalized  wave  heights  are  defined  by  X  =  H/Hq, 

X.^  =  H^/HQ  and  X2  =  H2/H0-  The  artifice  of  spreading  breakers 
over  a  range  partly  represents  the  inherent  variability  of 
breaker  heights,  and  partly  compensates  for  the  simplifica¬ 
tion  of  using  a  single  wave  period  in  the  estimation  of 
breaker  height.  Broken  waves  generally  have  some  height 
smaller  than  X^.  Since  no  theory  is  available  for  describing 
waves  after  they  have  broken,  the  heights  of  broken  waves 
are  assumed  to  be  redistributed  across  the  range  and  to  be 
proportional  to  the  unbroken  waves  as  was  done  by  Kuo  and 
Kuo  (1975)  .  Therefore,  the  conditional  probability  density 
function  for  all  heights  is  calculated  by 

P_(X) 

p  (X)  =  -  (16) 

X1 

/  p  (X) dx 
0  r 
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(13) 

(13) 

(15) 


where : 


1  -  -a  X, 

/  p  (X) dX  =  1-  { 1  +  at  X  (X.  -X,)}e  X,  (17) 

0  .  x  * 

is  a  constant  of  proportionality  applied  to  Pr(X)  to  normal¬ 
ize  the  pdf.  In  equation  (17),  the  constant,  a,  is  equal  to 

/7/ks. 

The  breaker  height  is  estimated  using  the  following 
formula,  which  is  an  approximate  expression  for  Goda's 
breaker  index  (1970)  based  on  laboratory  data, 

H  L»  H 

=  A  =2(1  -exp  [-1.5  Al  +K  tanp8)  )}  ,  (18) 

“o  o  o  ^o 

where  tan  6  denotes  the  bottom  slope  and  the  coefficients 
are  assigned  the  following  values  for  best-fitting  to  the 
index  curves, 

A  =  0.17,  K  =  15,  and  p  =  4/3. 

The  range  of  breaker  height,  X.^  -X2,  * s  calculated 
by  assigning  the  following  values  for  A: 

A^  =  0.18  for  X^ , 

A2  =  J  A1  =  0,12  for  X2* 

The  upper  limit  of  A^  was  selected  by  considering  the  varia¬ 
bility  of  breaker  heights,  whereas  the  lower  limit  of  A2  was 
chosen  simply  as  two-thirds  of  A^.  The  coefficients  used 
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are  based  on  matching  laboratory  data  taken  on  a  1/10  and 
1/50  beach  slope  and  several  field  experiments.  The  Goda 
model  is  applied  here  to  the  experimental  data  described 
below. 

5.  Summary 

The  common  idea  of  these  studies  is  to  cut-off  the 
portion  of  wave  height  distribution  beyond  the  breaker 
height,  which  is  controlled  by  the  water  depth  and  other 
factors.  The  methods  differ  in  the  techniques  of  cut-off 
and  the  formulae  used  to  define  breaker  heights. 
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III.  EXPERIMENT 


Experiments  were  conducted  at  Torrey  Pines  Beach,  San 
Diego,  California,  during  November  1978,  as  part  of  the 
Nearshore  Sediment  Transport  Study.  At  this  site  there  is 
a  gentle  sloping,  moderately  sorted,  fine-grained  sandy  beach. 
The  beach  profile  shows  no  well-developed  bar  structure  and 
is  remarkably  free  from  longshore  topographic  inhomogenei¬ 
ties.  Winds  during  the  experiments  were  light,  and  variable 
in  direction.  Shadowing  by  offshore  islands  and  offshore 
refraction,  limits  the  angles  of  wave  incidence  in  10m  depth 
to  less  than  15°.  During  the  experiments,  significant  off¬ 
shore  wave  heights  varied  between  60  and  160  cm.  The  condi¬ 
tion  of  nearly  normally  incident,  spilling  (or  mixed  plung¬ 
ing-spilling)  waves,  breaking  in  a  continuous  way  across  the 
surf  zone,  prevailed  during  most  of  the  experiments. 

A.  INSTRUMENTS 

An  extensive  array  of  instruments  was  deployed  to  study 
nearshore  wave  dynamics.  Measurements  described  here  are 
from  sensors  located  on  an  offshore  transect  from  10  m  depth 
to  across  the  surf  zone  (Fig.  1) .  The  sensors  were  of  three 
types:  pressure  (P) ,  current  (C)  and  surface-piercing-staff 

(W)  . 

The  pressure  sensors  were  Stathem  temperature -compensated 

2 

transducers  with  dynamic  range  of  either  912-2316  g/cm  or 
2 

912-3720  g/cm  .  They  were  statically  precalibrated  and 


postcalibrated  by  being  lowered  into  a  salt-water  tank  and 
were  found  quite  linear;  the  gains  differed  by  less  than  2% 
between  calibrations. 

Current  meters  were  two-axis ,  Marsh-McBirney  electro¬ 
magnetic,  spherical  (4  cm  diameter)  probes,  with  a  three- 
pole  output  filter  at  4  Hz.  Precalibration  and  postcalibra¬ 
tion  of  current  meters  showed  little  change  in  replicate 
runs  with  steady  or  oscillating  velocity  fields.  The  uncer¬ 
tainty  associated  with  using  a  single  gain  factor  for  all 
frequencies  is  roughly  estimated  at  ±5%  in  amplitudes  (10% 
in  variances) . 

The  wave  staffs  were  dual  resistance  wires  with  low 
noise,  high  resolution,  and  good  electronic  stability.  The 
accuracy  of  the  wave  staffs  was  about  ±3%  based  on  repeata¬ 
bility  of  gain  calibrations  measured  in  the  laboratory  and 
in  situ. 

B.  DATA  ANALYSIS 

Sea  surface  elevation  and  wave  velocity  components  were 
retrieved  from  sensors  by  telemetering  to  shore  and  there 
recorded  on  a  special  receiver/tape  recorder,  described  in 
detail  by  Lowe  et  al . ,  (1972)  .  The  sampling  rate  was  64 
samples/s  which  was  reduced  to  2  samples/s  by  digital  low- 
pass  filtering.  Record  lengths  of  approximately  68  minutes 
from  each  data  set  were  analyzed. 

It  was  desired  to  examine  only  the  sea-swell  band  of 
frequencies  between  0.05  to  1.0  Hz  (20  to  1  s  periods).  The 
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data  were  first  linearly  detrended  to  exclude  effects  of  the 
rising  and  falling  tides.  The  da'.a  were  then  high  pass 
filtered  with  a  cut-off  frequency  of  0.05  Hz  (20  s  period). 

The  high  pass  filter  used  a  Fast  Fourier  Transform  algorithm 
to  obtain  the  amplitude  spectrum  of  the  entire  68  minute 
record.  The  Fourier  coefficients  corresponding  to  0  to 
0.05  Hz  were  used  to  synthesize  a  low  frequency  time  series 
which  was  subtracted  from  the  wave  record.  The  limiting  high 
frequency  (Nyquist)  was  1.0  Hz. 

The  energy  density  spectra  were  calculated  in  a  similar 
manner  for  wave  and  velocity  measurements  using  the  Fast 
Fourier  Transform  (FFT)  algorithm.  A  cosine-squared  taper 
data  window  was  applied  to  the  time  series  to  minimize  leakage. 

The  highest  maximum  and  lowest  minimum  of  the  surface 
elevation  within  a  period  interval  defined,  respectively, 
the  crest  and  trough  of  a  wave.  A  wave  height  H  is  defined 
as  the  total  range  of  £(t)  in  that  interval,  the  time  between 
two  consecutive  zero-up-crossings  of  c(t)  (see  Fig.  2). 

Since  the  average  wave  period  was  about  14  sec,  the  total 
number  of  waves  in  the  68  min.  record  was  about  300,  which 
gives  reasonable  wave  statistics.  The  height  statistics  of 
mean  wave  height  H,  root  mean  square  wave  height  H  , 
significant  wave  height  H (average  of  the  heights  of  the 
1/3  highest  waves) ,  and  H^^q  (average  of  the  heights  of  the 
1/10  highest  waves) ,  are  calculated  from  the  ordered  set  of 
wave  heights. 
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The  pdf  and  cumulative  distributions  of  wave  heights 
were  calculated.  The  heights  were  normalized  using  the 


deep  water  root  mean  square  value.  Theoretical  probability 
distributions  were  calculated  using  the  Goda  model  and  com¬ 
pared  with  measured  distributions. 

A  deep  water  reference  wave  height  was  calculated  by 
measuring  the  energy  using  current  meter  C9  located  at  about 
4  m  depth  and  backing  the  energy  out  to  deep  water.  Kinetic 
energy  spectra  were  calculated  from  the  measured  horizontal 
velocity  spectra,  Su ( f )  and  Sv(f);  linear  theory  transfer 
functions  were  used  to  integrate  the  spectra  over  the  water 
column,  so  that  the  average  kinetic  energy  is 

<KE  (f )  >  =  |HRE(f)  |2(Su(f)  +Sv(f)  ]  (19) 

2 

where  |HRE(f) |  is  the  transfer  function  that  relates  the 
velocity  spectrum  components  to  the  kinetic  energy. 


|HKE(f) 


sinh  2kh 


cosh  k(h  +  z  ) 
m 


where  z  is  the  measurement  elevation.  Guza  and  Thornton 
m 

(1980)  showed,  for  these  same  experiments,  that  using  linear 
theory  transfer  functions  to  calculate  average  kinetic  energy 
gave  reasonable  results.  To  first  order  in  energy,  the 
average  potential  energy  equals  the  average  kinetic  energy. 
The  potential  energy  in  deep  water  is  obtained  by  applying 
linear  shoaling  transformation 


F 


<PE (f ) > 


deep  water 


IVf)  I  <pE(f)>4m  » 


(21) 


where  H  (f)  =  VC /C  ,  is  the  linear  shoaling  coefficient. 
5  go  g 

The  rms  wave  height  is  related  to  the  PE  and  KE  by 


<PE>  =  ~  p  g  H*  =  <KE>  =  /  <KE  (f )  >df  ,  (22) 

o  rms  y g 

from  which  the  deep  water  rms  wave  height,  denoted  hence¬ 
forth  by  Hq/  can  be  found. 

To  take  advantage  of  the  large  number  of  current  meters, 
current  data  were  used  to  infer  wave  heights .  The  velocity 
signals  were  convolved  using  linear  wave  theory  to  obtain 
surface  elevations.  The  complex  Fourier  spectra  of  the 
horizontal  velocity  components  U(f) ,  V(f)  were  first  calcu¬ 
lated  and  vectorially  added.  The  complex  surface  elevation 
spectrum,  X(f),  was  calculated  applying  the  linear  wave 
theory  transfer  function,  H(f) 


X(f)  =  H(f)-V(f) 


(23) 


where 


H  (f ) 


sinh  kh _ 

u>  cosh  k  (h  +  z  ) 
m 


(24) 


The  complex  surface  elevation  spectrum  was  then  inverse 
transformed  to  obtain  the  surface  elevation  time  series  from 


which  the  wave  height  distribution  is  calculated.  The  entire 
68  min.  record  was  convolved  at  one  time  in  order  to  minimize 
the  end  effects  which  result  in  spectral  leakage. 


IV.  RESULTS 


A.  TYPICAL  SPECTRA 

A  broad  range  of  wave  and  weather  conditions  were  en¬ 
countered  during  the  experiments.  Typical  velocity  spectra 
for  two  days  (Fig.  3,  lower  panel)  include  an  example  of 
very  narrow  band  spectra  calculated  for  November  20  for  the 
current  meters  C22x  and  C23x  which  straddled  the  mean  breaker 
line.  In  shallow  water  depths,  the  waves  generally  become 
more  "peaky",  resulting  in  increased  energy  at  the  harmonics. 
The  presence  of  strong  harmonics  in  the  spectra  indicates 
the  importance  of  nonlinearities  of  the  waves  in  shallow 
water.  The  spectral  energy  level  decreases  at  all  frequen¬ 
cies  except  at  the  very  lowest  from  the  deeper  instrument 
C22x  (heavy  line)  to  the  shallower  instrument  C23x  due  to 
breaking.  The  other  typical  spectra  (Fig.  3,  upper  panel) 
are  an  example  of  combined  sea  and  swell  with  a  narrow  band 
of  energy  at  swell  frequencies,  but  with  broad  band  energy 
at  higher  frequencies. 

B .  HEIGHT  STATISTICS 

The  wave  height  statistics  for  six  days,  inferred  from 

wave  staff  and  current  meter  data,  are  presented  in  Table  II. 

The  statistical  parameters  listed  are:  root-mean-square 

heights,  Hrms/  significant  wave  heights,  H^y3,  average  of 

the  heights  of  the  1/10  highest  waves,  H^^g,  and  the  maximum 

height,  H  .  These  parameters  were  obtained  from  68  min. 
max 
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records.  The  reference  wave  height  HQ,  obtained  by  backing 
the  energy  measured  at  current  meter  C9  out  to  deep  water, 
the  frequency  peak  of  the  spectrum  and  the  depth  for  each 
instrument,  are  also  presented.  Most  of  the  wave  staff 
measurements  were  made  inside  the  surf  zone,  after  the  waves 
have  started  breaking. 

Shoaling  effects  are  observed  in  the  data  of  Table  II. 

The  H^g,  and  ’  i-ncrease  as  depth  decreases, 

until  they  reach  the  breaking  point  after  which  the  wave 
heights  decrease  as  the  depth  decreases,  due  to  breaking. 

The  range  of  depths  of  the  instruments  is  from  about 
570  cm  to  40  cm.  The  average  peak  frequency  (peak  frequency) 
varied  little  during  the  experiments  and  was  about  0.07  sec  ^ 
(Table  II) .  The  relative  depth  for  the  waves  at  the  peak 
frequency  is  h/L  <  1/25  in  all  cases  so  that  they  can  be 
considered  shallow  water  waves. 

C.  COMPARISON  OF  EMPIRICAL  WITH  MODEL  DISTRIBUTIONS 

The  Goda  model  is  compared  with  measured  data.  The  model 

is  first  run  using  Goda's  original  coefficients.  This  model 

results  in  over-prediction  of  the  H^g.  The  coefficients 

are  changed  to  optimize  the  model's  description  of  the  wave 

height  distribution  qualitatively,  and  H^g  quantitatively. 

The  H  parameter  was  calculated  from  the  second  moment  of 
rms 

the  wave  height  distribution  and  is,  therefore,  a  more  sensi¬ 
tive  parameter  to  describe  the  shape  of  the  distribution 
(particularly  the  tail)  than,  say,  first  moments  such  as  the 
mean  H,  ^1/3*  ^1/10* 
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A  comparison  (Fig.  4)  was  made  between  the  empirical 
wave  height  distribution  normalized  with  the  reference  deep 
water  wave  height  and  Goda 1 s  model  distribution  for  November 
20.  Starting  in  deep  water,  it  is  observed  that  the  model 
fits  quite  well  the  empirical  distribution  obtained  from  the 
pressure  sensor  P7.  it  is  clearly  seen  that  the  wave  heights 
in  deep  water  are  essentially  Rayleigh  distributed.  In 
shallow  water  the  model  over-predicts  the  tails  and  under¬ 
predicts  the  peaks.  The  smaller  the  depth  the  greater  the 
errors  (Figs.  4  and  5).  Due  to  this  over-prediction  at  the 
tails  of  the  distribution  the  rms  wave  heights  obtained  by 
the  model  are  larger  than  the  measured  ones. 

To  obtain  a  better  agreement  between  the  measured  and 

the  calculated  rms  wave  heights,  the  higher  limit  of  breaking 

and  the  range  of  breaking  of  the  model  are  changed.  Goda 

(1975)  calculated  the  shoaling  using  the  nonlinear  theory 

of  Shuto  (1974)  (see  equations  10,  11  and  12) .  For  all  the 

ranges  of  wave  heights ,  frequencies  and  depths  used  here , 

2  2 

the  values  of  gHT  /h  are  calculated  and  in  all  cases  analyzed, 
fall  in  the  third  category,  equation (12)  .  This  law  was  used 
to  calculate  the  nonlinear  shoaling  coefficient  Kg .  Goda 
assumed  a  range  of  breaking  between  HI  and  H2  that  takes 
into  account  the  variability  of  breaker  heights  and  the  use 
of  a  single  frequency.  The  variable  breaker  height  is  a 
function  of  the  frequency,  depth  and  bottom  slope  (see  equation 
18) .  The  values  of  the  coefficients  used  by  Goda  in  his 
breaking  criterion  are  the  following:  Al  =  0.18,  A2  =  2/3  Al, 
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K  *  15  and  p  =  4/3,  empirically  assigned  to  give  the  best- 
fit  with  observed  wave  heights.  The  coefficients  used  here 
to  give  a  qualitatively  better  fit  with  the  measured  dis¬ 
tributions  and  quantitatively  fit  with  are:  A1  =  0.136, 

A2  »  1/2  Al,  K  =  20;  p  was  left  equal  to  the  previous  value. 
The  most  sensitive  coefficient  is  Al  which  was  determined 
first.  Values  of  Al  were  calculated  for  various  distributions 
by  first  obtaining  the  higher  breaker  limit  (XI)  from  the 
measured  empirical  distributions  as  indicated  by  the  maximum 
value  of  the  distribution,  and  then  calculating  Al  using  equa¬ 
tion  18 .  The  values  of  Al  were  then  averaged  to  give  a 
single  representative  value  for  all  distributions.  The 
coefficient  A2  was  chosen  as  1/2A1,  which  results  in  a  more 
symmetrical  distribution  as  indicated  by  the  results.  The 
coefficient  K,  which  weights  the  slope  of  the  beach,  was 
determined  by  trial-and-error  testing  of  Goda's  model  for  a 
variety  of  values  looking  for  the  best  fit  of  all  the  distri¬ 
butions.  The  model  is  not  very  sensitive  to  changes  in  K. 

The  comparison  (Fig.  6)  of  the  empirical  distribution 
in  deep  water  for  November  20  with  the  distribution  obtained 
applying  the  model  used  the  new  coefficients.  In  deep  water, 
there  is  no  notable  difference  from  the  original  model's 
results  (Fig.  4) .  In  shallow  water,  the  predicted  values 
obtained  with  the  modified  model  fit  much  better  the 
empirical  distributions  than  those  obtained  applying  the 
original  Gcda  model. 
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As  stated  earlier,  the  coefficients  Goda  originally  speci¬ 
fied  were  based  on  matching  laboratory  data  for  a  two  beach 
slopes  of  1/10  and  1/50.  The  model  was  then  applied  to  some 
field  data,  but  unfortunately  the  beach  geometry  (slopes) 
were  not  given.  The  Torrey  Pines  Beach  is  approximately 
1/50  at  the  beach  face  and  the  wave  climate  is  characterized 
by  long  period  (~14  sec)  swell.  The  reason  for  differences 
in  the  model  coefficients  needed  to  fit  this  data  set  com¬ 
pared  with  Goda's  suggested  coefficients  is  not  known. 

D.  COMPARISON  OF  RMS  WAVE  HEIGHTS 

The  model  calculated  H  values  and  measured  values  of 

rms 

H _  calculated  directly  from  the  wave  heights  are  used  to 

rms  3 

test  how  well  the  model  works.  Table  III  shows  the  measured, 

Rayleigh,  and  the  calculated  Goda  and  modified  Goda  root- 

mean-square  wave  heights  corresponding  to  the  wave  staffs 

and  current  meters  for  six  different  days. 

Goda's  model-predicted  H _  values  are,  in  general, 

rms 

larger  than  the  measured  ones  (Fig.  9) .  In  an  attempt  to 

explain  the  differences,  H  „  values  obtained  from  the 

rms 

model  were  plotted  against  depth,  deep  water  wave  height 

and  Ursell  number;  but  no  meaningful  correlation  was  obtained 

with  any  of  these  variables . 

Measured  H  and  calculated  values  assuming  the  wave 
rms 

heights  are  Rayleigh  distributed  so  that  Hrms  =  /8mo  were 
compared  (Fig.  8a) .  Measured  and  assumed  Rayleigh  signifi¬ 
cant  wave  heights  were  also  compared  (Fig.  8b)  .  The  comparisons 
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show  that  the  agreement  between  measured  and  Rayleigh  sta¬ 
tistics  are  good  for  both  small  and  large  wave  heights, 
inside  and  outside  the  surf  zone.  The  good  comparisons  with 
Rayleigh  statistics  suggests  that  the  wave  heights,  although 
decreased  by  breaking,  are  still  more  nearly  Rayleigh-dis¬ 
tributed  than  cut-off  Rayleigh,  as  suggested  by  Goda. 

The  change  of  the  rms  wave  heights  with  depth,  between 
the  measured  and  the  calculated  ones  obtained  by  the 

model  first  applying  nonlinear  (heavy  line)  and  then  linear 
shoaling  (light  line)  were  compared  (Fig.  10).  The  model  with 
nonlinear  shoaling  clearly  over-predicts  the  values,  while 
the  model  with  linear  shoaling  gives  reaonable  values,  com¬ 
pared  with  measurements.  The  nonlinear  shoaling  "blows  up" 
in  very  shallow  water  (<30  cm)  and  should  be  ignored. 

Measured  and  the  calculated  rms  wave  heights  with  depth 

obtained  by  applying  linear  and  nonlinear  shoaling  using  the 

modified  coefficients  were  also  compared  (Fig.  10,  lower 

panel) .  This  figure  illustrates  that  the  modified  Goda  model 

fits  better  the  data  than  the  original  does.  The  majority 

of  the  measured  H _  fall  between  the  two  curves  obtained  with 

rms 

nonlinear  and  linear  shoaling.  Based  on  the  choice  of  coeffi¬ 
cients,  applying  linear  shoaling  to  Goda's  model  can  give 
as  good,  or  better,  results  as  applying  nonlinear  shoaling. 

The  rms  wave  height  values,  obtained  using  the  modified 
model  coefficients,  were  plotted  against  depth,  deep  water 
wave  height  and  Ursell  number;  no  obvious  correlation  could 
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be  noted.  It  is  found  that  the  new  predicted  H  values 

rms 

compared  with  measured  H  have  an  error  of  less  than  ±20% 

rms 

at  all  depths  (with  the  exception  of  one  anomalous  point) 

(Fig.  11) . 

E.  WAVE  HEIGHT  DISTRIBUTIONS  USING  CURRENT  METERS 

As  described  earlier,  the  surface  elevations  were  derived 
by  linearly  convolving  the  velocity  records  and  wave  height 
distributions  calculated.  The  basis  for  applying  this  analy¬ 
sis  is  the  earlier  work  of  Guza  and  Thornton  (1980)  where 
they  showed,  for  this  same  data  set,  that  linear  theory  spec¬ 
tral  transformations  could  be  used  to  calculate  surface  ele¬ 
vation  standard  deviations  either  from  pressure  meters  or 
current  meters  with  less  than  a  20%  error,  and  typically  less 
than  10%.  Examples  of  the  derived  wave  height  distributions 
are  considered  for  November  17  (Fig.  12).  These  measurements 
were  made  just  outside  the  surf  zone,  at  about  the  breaker 
point  and  inside  the  surf  zone  (current  meters  C23,  C37  and 
C40  respectively).  In  general,  the  model  overestimates  the 
velocity  derived  wave  height  distributions  more  than  the 
direct  measurements .  The  reason  for  the  discrepancy  is  that 
linear  wave  theory  underestimates  the  surface  elevations  in 
convolving  the  velocities,  particularly  in  the  crest  region 
of  the  waves.  In  other  words,  linear  theory  does  not  account 
for  the  finite  amplitude  of  these  highly  nonlinear  waves. 


a 
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F.  COMPARISON  OF  MODEL  AND  MEASURED  CUMULATIVE  DISTRIBUTIONS 
The  cumulative  exceedance  of  wave  height  distributions 
normalized  with  rms  deep  water  wave  height  were  calculated 
applying  the  modified  Goda  model  (using  nonlinear  shoaling) 
and  plotted  with  the  measured  cumulative  distribution  for 
comparison.  The  cumulative  exceedance  distribution  empha¬ 
sizes  information  in  the  tail  of  the  distribution.  For  an 
example  in  shallow  water,  inside  the  surf  zone  (116  cm  depth) , 
there  is  a  good  agreement  between  the  measured  and  predicted 
distributions  with  a  slight  underprediction  by  the  model  in 
the  tail  (Fig.  13) .  With  sensors  located  just  outside  the 
surf"  zone,  under-prediction  of  the  tail  is  larger  and  over¬ 
prediction  in  the  middle  range  occurs  (Fig.  14) .  In  general, 
there  is  a  better  agreement  between  the  two  distributions  well 
inside  the  surf  zone,  e.g.,  wave  staffs  W38  and  W41,  than 
those  (W21  and  W29)  which  were  generally  either  at  breaking 
or  just  outside  the  surf  zone,  the  most  nonlinear  wave 
region. 
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V.  CONCLUSIONS 


It  is  confirmed  that  the  wave  heights  in  deeper  water 
(7  m)  are  Rayleigh  distributed. 

Goda's  model,  using  the  empirical  coefficients  originally 
suggested  on  the  basis  of  laboratory  and  poorly  specified 
field  measurements,  over-predicts  the  tail  of  the  distribution 
and  under-predicts  the  peaks.  As  a  consequence,  the  predicted 
rms  wave  heights  are  larger  than  the  measured.  As  the  depth 
decreases,  the  errors  of  the  predicted  distributions  increase. 
To  get  a  better  fit  with  the  measured  wave  height  distribu¬ 
tions  ,  the  coefficients  in  the  breaking  criterion  used  in 
Goda's  model  were  modified.  The  values  for  Goda's  breaking 
criterion  giving  the  best  fit  to  the  measured  data  of  this 
study  are: 

Al  =  0.136,  A2  =  1/2  A1  and  K  =  20  . 

The  percentage  of  error  between  the  measured  and  the  pre¬ 
dicted  H _  values  from  the  model  with  these  coefficients  is 

rms 

les  than  ±20%. 

Linear  shoaling  was  found  to  be  as  good  as  nonlinear 
shoaling  in  applying  Goda's  model  across  the  surf  zone. 

Good  comparisons  were  obtained  between  empirical  and 
Rayleigh-derived  statistics.  This  indicates  that  the  wave 
heights,  although  decreased  by  breaking,  are  still  more 
Rayleigh-distributed  than  the  cut-off  Rayleigh-distributed 
as  suggested  by  Goda. 
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Goda's  model  was  tested  here  with  a  large  amount  of 
field  data  for  a  variety  of  wave  conditions  on  a  1/50  beach 
slope.  Further  comparisons  should  be  made  for  a  variety  of 
beach  slopes  and  wave  climates  in  order  to  test  the  general 
applicability  of  the  model. 


TABLE  II.  WAVE  HEIGHT  STATISTICS 


(all 

values 

in  cm) 

Date 

Inst 

h 

H 

o 

f 

H 

rms 

Hl/3 

Hl/10 

H 

max 

Nov 

4 

W41 

82 

42.5 

.0703 

33.2 

44.4 

53.2 

74.3 

W38 

125 

42.5 

.0703 

44.8 

61.5 

73.4 

98.1 

W21 

177 

42.5 

.0703 

50.6 

71.8 

90.9 

131.3 

W29 

225 

42.5 

.0703 

56.7 

81.3 

113.8 

188.3 

Nov 

10 

W21 

170 

66.1 

.  0632 

61.0 

85.5 

107.4 

140.4 

Nov 

17 

W41 

92 

44.8 

.0729 

37.8 

50.7 

59.4 

79.5 

W38 

141 

44.8 

.0729 

52.0 

72.1 

86.9 

111.5 

W21 

197 

44.8 

.0729 

54.1 

76.2 

94.1 

139.2 

W29 

209 

44.8 

.0729 

52.0 

72.8 

93.7 

154.0 

Nov 

20 

W38 

116 

52.4 

.0666 

35.9 

49.3 

57.9 

83.8 

W21 

153 

52.4 

.0666 

48.6 

68.1 

84.2 

120.3 

W29 

182 

52.4 

.0666 

73.1 

110.2 

143.0 

202.7 

Nov 

20 

C42 

39 

52.4 

.0703 

18.0 

27.0 

33.8 

49.1 

C39 

102 

52.4 

.0703 

35.6 

50.5 

59.8 

82.6 

C37 

116 

52.4 

.0703 

36.5 

53.6 

62.1 

81.8 

C36 

142 

52.4 

.0703 

40.1 

60.4 

69.8 

89.5 

C23 

147 

52.4 

.0703 

51.6 

80.2 

98.0 

118.9 

C22 

188 

52.4 

.0703 

60.8 

96.1 

127.7 

159.3 

C19 

250 

52.4 

.0703 

68.9 

105.1 

149.0 

210.4 

C15 

355 

52.4 

.0703 

58.7 

88.5 

125.5 

217.9 

C09 

571 

52.4 

.0703 

55.8 

84.0 

115.2 

192.4 

Nov 

24 

W38 

65 

36.4 

.0639 

10.1 

14.5 

18.4 

26.3 

W21 

86 

36.4 

.0715 

26.7 

36.6 

43.5 

55.8 

Nov 

18 

W41 

84 

55.4 

.0757 

33.8 

45.8 

55.4 

109.1 

W21 

195 

55.4 

.1552 

61.4 

87.2 

103.9 

131.7 

W29 

198 

55.4 

.0756 

63.9 

91.3 

117.2 

166.3 

TABLE  III. 


MEASURED  AND  CALCULATED  RMS  WAVE  HEIGHTS 
OBTAINED  WITH  GODA'S  MODEL  AND  MODIFIED 
GODA'S  MODEL  (all  values  in  cm) 


RMS 

;  WAVE 

HEIGHTS 

SIG. 

HEIGHTS 

Date 

Inst. 

Meas . 

Ray. 

Goda 

Mdf .G. 

Meas. 

Ray 

Nov 

4 

W41 

33.2 

33.7 

34.6 

30.9 

44.4 

46.9 

W38 

44.8 

46.2 

55.9 

47.5 

61.5 

63.4 

W21 

50.6 

52.3 

67.3 

51.1 

71.8 

71.6 

W29 

56.7 

58.2 

65.3 

58.5 

81.3 

80.2 

Nov 

10 

W21 

61.0 

64.9 

85.3 

64.5 

85.5 

86.3 

Nov 

17 

W41 

37.8 

37.9 

48.6 

34.7 

50.7 

53.5 

W38 

52.0 

52.2 

66.1 

53.4 

72.1 

73.5 

W21 

54.1 

55.6 

70.3 

56.0 

76.2 

76.5 

W29 

52.0 

53.8 

69.8 

58.1 

72.8 

73.5 

Nov 

W38 

35.9 

36.2 

50.3 

43.9 

49.3 

50.8 

W21 

48.6 

50.3 

68.3 

44.9 

68.1 

68.7 

W29 

73.1 

67.5 

77.6 

53.3 

110.2 

103.4 

Nov 

20 

C42 

18.0 

18.2 

21.4 

18.0 

27.0 

25.4 

C39 

35.6 

36.1 

54.0 

38.5 

50.5 

50.3 

C37 

36.5 

36.3 

60.3 

43.9 

53.6 

51.6 

C36 

40.1 

40.0 

70.1 

53.8 

60.4 

56.7 

C23 

51.6 

50.0 

71.5 

55.7 

80.2 

73.0 

C22 

60.8 

57.5 

79.1 

54.9 

96.1 

86.0 

C19 

68.9 

65.0 

78.3 

67.5 

105.1 

97.4 

CIS 

58.7 

60.9 

72.0 

64.0 

88.5 

83.0 

C09 

55.8 

60.7 

64.0 

55.6 

84.0 

79.0 

Nov 

24 

W38 

10.1 

10.7 

26.7 

24.4 

14.5 

14.3 

W21 

26.7 

27.5 

44.6 

32.7 

36.6 

37.8 

Nov 

18 

W41 

33.8 

35.5 

33.4 

31.4 

45.8 

47.8 

W21 

61.4 

59.7 

59.6 

50.8 

87.2 

86.8 

W29 

63.9 

61.7 

80.7 

57.5 

91.3 

90.4 
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where 


X  is  a  crest 

*  is  a  trough 

0  is  the  zero-up  crossing 


Figure  2.  Definition  sketch  of  zero-up-crossing 
wave  heights . 
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‘  LOG  SPECTRAL  DENSITY  CM*m2/SEC 

LOG  SPECTRAL  DENSITY  CM«m2/SEC 

J.O  2.0  3.0  4.0  5.0  o1.0  2.0  3.0  4.0  5.0 


******* 


NOV  4  FILE  1 


FREQUENCY  (HZ) 

Figure  3.  Typical  spectra  measured  during  the  experiments. 
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Empirical  distribution  of  wave  heights  compared 
with  those  predicted  with  Goda's  model,  starting 
in  deep  water  (P7)  and  going  into  shallow  water 
(W21,  W38) ,  20  November  1978. 
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Figure  5. 


Empirical  distribution  of  wave  heights  compared 
with  those  predicted  with  Goda's  model,  starting 
deep  water  (W21)  and  going  into  shallow  water 
38,  W41) ,  17  November  1978. 
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Figure  6.  Empiric 
with  th 
model. 
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Figure  8a.  Range  of  measured  and  Rayleigh  root-mean- 
square  wave  heights. 
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Figure  8b.  Range  of  measured  and  Rayleigh  significant 
wave  heights . 
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Figure  12.  Empirical  distributions  of  wave  heights 

obtained  from  current  meters  (C23,  C37  and 
C40)  compared  with  predicted  wave  heights 
calculated  with  the  modified  model, 

17  November  1978. 
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Figure  13.  Comparison  of  measured  cumulative  exceedance 
distributions  with  predicted  distributions 
(modified  model) ,  20  November  (W38)  and 
17  November  1978  (W41) . 
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Figure  14.  Comparison  of  measured  cumulative  exceedance 
distributions  with  predicted  distributions 
(modified  model) ,  20  November  (W21)  and 
17  November  1978  (W38) . 
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